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Abstract 

A superconducting nanoscopic cylinder, with radius R — 4.0^ and height D — 4.0^ is submitted 
to an applied field along the cylinder axis . The Ginzburg-Landau theory is solved in three- 
dimensions using the simulated annealing technique to minimize the free energy functional. We 
obtain different vortex patterns, some of which are giant vortices and up to twelve vortices are able 
to fit inside the cylinder. 
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Recently, a diode based on vortex physics was proposed'lj. It exhibits new features 
and opens new prospects for flux hne circuit implementations. Potential candidates for 
resistor, capacitor and inductor based on vortex physics were also considered and they all 
rely on the choice of sample geometry. Vortex patterns have been previously studied in 
several geometries with nanoscopic dimensions, such as disks [2], triangles [3] and squaresP]. 
They were theoretically studied using the Ginzburg-Landau (G-L) theory but only in the 
two-dimensional limit, where the sample thickness is smaller than ^. This allowed them to 
assume that the order parameter is constant along the third direction deflned by the fleld 
reducing the dimensionality of the problem. The other limit of an inflnite long cylinder was 
studied in [5]. Here we consider the intermediate case of a nanoscopic cylinder with radius 
R = 4.0(^ and height D = 4.0^ using the full three-dimensional approach and obtain the 
vortex states as a function of the applied fleld perpendicular to its surface. 

We minimize numerically the G-L free energy functional using the simulated annealing 
procedure, discussed in Ref. [6] and in references quoted therein. Every vortex state is 
characterized by an integer number which is the total vorticity inside the cylinder. In some 
magnetic fleld range we flnd for increasing fleld and constant vorticity that the vortex evolves 
from a set of individual vortices with unitary vorticity to another set containing vortices 
with vorticity larger than one, the so called giant vortex. However for the same conditions 
of increasing fleld and constant vorticity the pattern becomes metastable as another vortex 
state becomes lower in free energy. This change of stability takes place at the matching 
fleld jzl- In flgure n] we show the free energy versus the applied fleld and in flgure [21 the 
corresponding magnetization. The lowest energy curve corresponds to the thermodynamic 
most stable pattern, the Meissner phase with no vortices, but only up to the flrst matching 
fleld. Above this fleld the entrance of a single vortex is favored, though it is still possible 
to keep the system in the Meissner phase but now as a metastable conflguration. The 
small difference in energy around the matching fleld is able to drive the state to a magnetic 
hysteresis. In flgure we flnd Multi Vortex States (MVS) and Giant Vortex States (GVS) 
for flelds below a threshold value H/Hc2 = 0.96 and only GVS above it. This result is 

n 

different from those in Ref [8] where they found only GVS as ground states (MVS appeared 
as metastable states). The red and blue lines in the magnetization, flgure [21 present the 
full sweep fleld up and down, respectively. The magnetization curve is truly a collection 
of several independent lines, each associated to a distinct vorticity which do not cross each 
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other. As the apphed field is swept the vortex pattern jumps from one curve to the next 
one causing the vorticity to jump with one unit and leading to hysteretic behavior. Figures 
2B and 3C show the Cooper pair density of two different states corresponding to the field 
H/Hc2 = 0.785. The state with L = 4 was obtained with increasing the field(red curve) 
and the state L = 6 with decreasing the field. We find twelve independent magnetization 
lines which implies that a total of twelve vortices fit inside the cylinder. This result is in 
fair agreement with the eleven vortices found by Baelus et al. for a disk with the same 
radius, thickness 0.1^ and = 0.28. Their study is based on a different method, the time 
evolution of the corresponding two-dimensional system. The maximum number of vortices 
inside a cylinder follows from the ratio between two areas, the disk surface and the vortex 
core: Ug = /^^ = 16. This value overestimates the true result found here because it does 
not consider the repulsion among vortices. The first vortex appear at the field H/ Hr2 — 0.21 
and the normal phase is reached at the applied field Hcs/Hc2 = 1.48. For a disk [8] larger 
values were found for these critical fields, namely, 0.38 and 1.95 in units of Hc2 (which is a 
consequence of the smaller expulsion of the magnetic field in a thin disk). 
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FIG. 1: Free energy of the GVS (solid curves) and the MVS (dashed curves) versus field for a 
cylinder of radius R = 4.0^ and thickness D = 4.0^. The labels inside the curves indicate the 
vorticity of the system. In the grey region only GVS were observed. In the white region both MVS 
and GVS were observed. 
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FIG. 2: (A) Magnetization versus applied field for the cylinder. The insets are the two dimensional 
|^^(r)p = cte plots for H/Hc2 = 0.785 with vorticity 4 and 6. (B) and (C) show the three dimension 
Cooper pair density isosurface. 



5 



